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We present an approximate gap equation for different crystalline structures of the LOFF phase of 
high density QCD at T=0. This equation is derived by using an effective condensate term obtained 
by averaging the inhomogeneous condensate over distances of the order of the crystal lattice size. The 
approximation is expected to work better far off any second order phase transition. As a function of 
the difference of the chemical potentials of the up and down quarks, 8fi, we get that the octahedron 
is energetically favored from 8[i — Ao/\/2 to 0.95Ao, where Ao is the gap for the homogeneous 
phase, while in the range 0.95Ao — 1.32Ao the face centered cube prevails. At <5/i = 1.32Ao a first 
order phase transition to the normal phase occurs. 



I. INTRODUCTION 

The existence of an inhomogeneous superconductive phase characterized by a periodic structure of the gap parameter 
was hypothesized long ago in two separate and contemporary papers by Larkin and Ovchinnikov Q and by Fulde and 
Ferrel |2| and is therefore named LOFF phase. These original papers studied superconductors in presence of a strong 
magnetic field coupled to the spins of the conduction electrons. The corresponding interaction hamiltonian gives rise 
to a separation of the Fermi surfaces of the pairing electrons since they have opposite spins. If the separation is very 
high the pairing is destroyed and there is a transition from the superconductive to the normal state. In Q and 0] it 
was shown that a new state could be formed, close to the transition line. This state has the feature of exhibiting an 
order parameter, the gap, which is not a constant, but has a space variation whose typical wavelength is of the order 
of the inverse of the difference in the Fermi energies of the pairing electrons. The space modulation of the gap arises 
because the electron pair has non zero total momentum. In the simplest case the modulation is that of a plane wave. 
In the general case it is given by the sum of several plane waves, which might lead to a crystalline structure. In any 
case the LOFF phase breaks translational and rotational symmetries. 

Inhomogeneous superconductivity in metals has been the object of intense experimental investigations especially 
in the last decade, though the evidence is still inconclusive (for a review see |3|). Quite recently it has been also 
realized that at moderate density the mass difference between the strange and the up and down quarks at the weak 
equilibrium leads to a difference in their Fermi momenta, which renders in principle the LOFF state possible in color 
superconductivity Q (for general reviews of color superconductivity see 0, IE IE Eil El)- The study of the 
crystalline color superconductive phase of QCD is theoretically important for the understanding of the QCD phase 
diagram. It may be also relevant to the astrophysics of compact stars. As a matter of fact whereas at extreme densities 
one expects the color-flavor-locking phase to occur, at lower densities, within the inner core of compact stars, there 
is a better chance that the crystalline phase be realized. In these conditions a difference in the chemical potentials 
of the up and down quark might be generated by beta-equilibrium of the quark matter. While difficult to be studied 
in the laboratory this phenomenon might be therefore relevant in astrophysics, for example by explaining the sudden 
variations (glitches)of the rotation period of the pulsars (for a review of color superconductivity in compact stars see 
e -g-E2> EJ)- I n E3 a Ginzburg-Landau (GL) analysis of several crystalline structures was performed for QCD at 
T = 0. This approximation holds for small values of the gap parameter and therefore this analysis studies a small 
region of the phase diagram, i.e. the corner where the difference of the chemical potentials of the pairing quarks d/j, is 
near to 5/i2, the value where a second order phase transition from the LOFF to the nor mal p hase is assumed to exist. 
For the one-plane- wave structure 5/J.2 is not far from the Clogston-Chandrasekhar [l5|, [fj| point <5/ii = Ao/V2 (Ao 
is the value of the gap for the superconductive homogeneous phase). This is approximately the point where a first 
order transition occurs between the homogeneous Bardeen-Cooper-Schricffcr (BCS) and the inhomogeneous LOFF 
phases. In other cases Sfj,2 can be larger and the phase transition can be first-order. 

An effective field theory for the crystalline color superconductive phase of QCD was developed in 01 ■ The papers 
[l8| and studied the properties of the Nambu-Goldstone Bosons (NGB) associated to the breaking of the translation 
symmetry, whose properties were analyzed in |20j |. where an approximation scheme valid at the GL point (A — * 0) was 
developed. In the present paper we come back to this problem with the aim to implement an approximation scheme 
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for the gap equation along similar lines. We will find that an effective gap equation can be indeed written down, but 
it is expected to work better for large A, for example close to the BCS/LOFF transition point 6 fix- We will work at 
T = and for high densities, since we wish to study the LOFF phase in the context of QCD, where the physically 
relevant case is presumably the zero temperature high density inner core of compact stars. We will not consider a 
complete analysis, which is of course an endless task, but we shall limit our analysis to four crystalline structures. In 
Section nil we discuss the one plane wave as the prototype case. Then we consider those cases that various analyses 
[l4| . |2l| indicate among the favored forms of the condensate: the two-antipodal plane wave, the octahedron and the 
face-centered-cube, to be discussed in Sections IIIII and llVl The numerical results of Section Hvl are that, as a function 
of the difference 5p between the quark chemical potentials, the octahedron is energetically favored from 5p = Ao/v2 
to 0.95A , while in the range 0.95A — 1.32A the face centered cube prevails. After this point the normal phase is 
restored. 



II. THE EFFECTIVE GAP EQUATION FOR THE FULDE-FERREL (ONE-PLANE WAVE) STATE 

We begin our study by a review of one-plane wave, Fulde-Ferrel (FF) state. Although this case can be solved 
exactly, it is useful to consider it here in order to fix the notations and introduce some definitions to be used later on. 
We shall consider Cooper pairing of the massless quarks up u and down d, with chemical potential p u , pd- We define 
M = (Mu + Md)/2 and Sp, = \p u - p d \/2 < iu 

The gap equation, as discussed in Ref. [3J, can be derived using the following effective interaction Lagrangian 

C^—g^y^^fX"^. (1) 

Here g is a coupling constant, A a are color matrices and a sum over flavors is understood. 

Following the same steps as in Q we get that, in the mean field approximation, Ci in Q is substituted by 

Ccond = -\e a ^{^ A(r) + c.c.) + (L -> R) - -A(r)A*(r) , (2) 
£ g 

where i, j are flavor indices and a, (3 are color indices. In the FF state the total momentum of the Cooper pair is 2q 
and the condensate has the space-dependence of a single plane wave J2j: 

A(r) = Ae 2iqr . (3) 

It follows that at zero temperature, the gap equation for the FF state can be written as follows 

gp f dv f s d£ 



1 = — - / — — / s (1 - 9(-E u ) - 0(-E d )) , (4) 

where |v| = 1 and 5 is the ultraviolet cutoff that we assume equal to /i/2; moreover we assume S 3> q ~ 5fi. The 
following equations 



E u4 — ±5(i =F q • v + Ve 2 + A 2 (5) 
provide the quasi-particle dispersion laws and p is the density of states that in QCD with two flavor quarks is 

P = — ■ (6) 

We observe that 

1 - 6(-E u ) - 6{-E d ) = 9{E u )B{E d ) (7) 
so that the integration in Eq.Q is over the pairing region (PR), defined by 

PR = {(£,v) | E u > Oand£ d > 0} . (8) 
More explicitly the pairing region (PR) is defined by the condition (q = q/q) 



"] f Jf 2 + A 2 ] 
■ <v-q< Min ll, z q + \ 



Max] -1, z q - ^ +A ) < v- q < Min { 1, z q + } {<)) 
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with 



z g = ^ • (10) 

q 



Eq. (0} can be written in a way that will be useful later: 



that after energy integration gives 



A^fg/^ ■ (12) 



47T ,/n + A 2^ 



Here A e // = A e //(v • q, £) is defined as 



(- A for (f , v) G Pi? 
A e// =A9(E u )d(E d ) = I (13) 

t elsewhere . 

To the interaction term J2J in the lagrangian one should add the Lagrangian for free quarks. We adopt for its 
description the formalism of the High Density Effective Theory (HDET) (see [22, H2, [3, EE US , an d, for reviews, 
[TUlH 113) and write it as follows 

£ o = [^+ iV ■ + ■ d% l J - + {L —> R) . (14) 

V 

Here the sum represents an average over velocities. The velocity dependent left-handed fields ip± = ip± v are the 
positive energy part in the decomposition 

V-(z) = / % er ^"" I4>v(x) + *v(*)] (15) 
while \? v is the negative energy part which has been integrated out. %p v is given by 

^v(x) = e^P+VC*) = / /^e^-P+^W (16) 

and therefore contains the residual momentum I, corresponding to the decomposition of the quark momentum p = 
Hv + I, with v<- L = (0,v) and t» = lv = £. We also introduce = (l,v) and V M = (1, — v). P± are projectors 
defined by 

P ± = i(l± 7o7 -v) . (17) 
In this formalism we get from P)l (C = ioi is the charge conjugation matrix) 

£cond = ~ ^ eX P{* r ' a ( V i> V j' q)} £ ij £ a/33V'J i ;ia( x ) C '0-v j ;j/3(x) 

V i: Vj 

-(L->R)+h.c - -A(r)A*(r) , (18) 
5 

where 

a(vt, vj, q) = 2q - /ijVi - /LtjVj . (19) 
This Lagrangian can be put, in momentum space, in the following form 

Uu=~ ^A e// ei J -e^ 3 < i «(W-v; J -/3M)-(£^-R)+li-c.--AA*, (20) 
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by an average procedure (see below) and then we can obtain the gap equation lllj l. 
following basis for the fermion fields: 



It is useful to introduce the 



A=0 



V2 



1,2), V 



+ ,31 



+ ,32 



1>l 



(21) 



where cr^ are the Pauli matrices for A = 1,2,3 and do = 1 (similar expressions hold for tp^) and then we obtain that 
the propagator for the velocity dependent fermionic fields is given by 



DabW 



1 



A 2 



V-£S AB -A ab 
-A AB V-£5 AB 



V-lV-t 

The matrix Aab is as follows: Aab — (A or B = 4 or 5), and, for A, B = 0, 3: 



(22) 



A AB = A 



eff 



( 1 

0-10 
0-10 

Vo o o-i 



(23) 



AB 



To obtain the Lagrangian (|20ll one can perform a weighted average of the Lagrangian H18|) over a region of the size 
of the lattice cell. The weight function gu{r) will be defined below. We note that in the gap equation the relevant 
momenta are small with respect to the gap which is of the order of q. Therefore we may assume that the velocity 
dependent fields are slowly varying over regions of the order of the lattice size. This means that in the average we 
can treat them as constant, and in conclusion the average is made only on the coefficient exp{ir ■ a}. Therefore what 
we are computing is 



1(a) = ^ exp{ir • a}g R (r) 



(24) 



where the bracket means average over the cell. 
It is possible to choose gn(r) in such a way that 



/(a) = 5% 



2g 



(25) 



where 



1 for Ixl < — 
? 3, , } ' ' 2R 



^(x) = 



(26) 



elsewhere . 



An approximate expression of gi?(r) is given by 



g R {T) = n 





'■*qr k - 


sin 


. R \ 



k=l 



(27) 



For R/n w 1, gn is different from zero only in a region of the size of the unit lattice cell; the resulting integral / 
obtained by l|24() and l|27|l is reported (for the one-dimensional case and for R = 2tt/3) in figQ 

Independently of the exact form of ffi?(r), we will assume that the average procedure gives as a result the brick- 
shaped function 5r defined in 12tj|) . Clearly this assumption is a crucial ingredient of our approximation. 

In the /i — > oo limit by taking q along the z— axis, we get for the components x and y of a 



nq 



i.e. approximately (for Sfi <C //) 



(vi + V 2 )x,y\ < 



Rjl 



(28) 



(29) 
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FIG. 1: I{x) as obtained by Equations an d ( m one dimension) for q = 40 MeV and i? = 27r/3. 



From this in the high density limit it follows 

vj = -v a + 0(l//z) ■ (30) 

We use the x and y components of the 6r function to get rid of the integration over V2 in eq. Ijl8(l . The two factors 
ir/R arising from the x and y components are absorbed into a wave function renormalization of the quark fields, both 
in the kinetic and in the gap terms. The z component gives a factor 5n[h(v ■ q)], where the function h was computed 
in pcj . Since here we work in the limit Sfj, <C fi, we neglect corrections of order Sfi/fi and we use the asymptotic 
expression [20j : 

fe(v4) = l-Ar- (31) 
v • q 

After this average we get that the coefficient of the bilinear term in Eq. I|19|l is 

A6 R [h(vq)]. (32) 
This factor coincides with A e //(v • q, £) if one chooses 

a= ^-yql . (33) 

2^/eT^\h(v ■ q)| 

This is consistent with our hypothesis, R/n m 1, only if A is not too small, meaning that we should be far from a 
second order phase transition. 

Let us finally notice that Eq. H12|) can be obtained from Eq. Ijllfl by substituting the previous expression for R 
with 

_ njfrt- vq| 

2|£ |-|Mv-q)| ■ 1 J 

In fact, observing that in any case A e ff is equal to or A, at the pole we get back the expression Q33f) . Then written 
as in i|34[l. i? and analogously A e // become functions of the velocity and the energy; therefore our average should be 
better taken in the momentum space. 

III. THE EFFECTIVE GAP EQUATION AT T = FOR GENERIC CRYSTALLINE STRUCTURES 

Let us now consider the general case of P plane waves 

p 

A(r) = Y, A m e 2M ^ r . (35) 

m— 1 

We shall consider only the case A m = A for any m = 1, ...P and 

q m = (?n m , (36) 
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with n m unit vectors. Generalizing the results of the previous equations we substitute in the Lagrangian l|20|l and in 
the propagator (|2*2)l A e jf (v • n,£ ) with 

p 

A B (v, £ ) = A -ff ( v ' n ™' £ °) ■ ( 3? ) 

m=l 

In this equation we have made explicit the dependence on the energy £q arising from the average. Correspondingly 
the gap equation is 

pa _ jP fdv fdH A E (v,£ ) 

' 2 J An J 2ir £q £jj A|(v,* ) 1 j 

which generalizes Eq. The origin of the factor P on the l.h.s. of this equation is as follows. The Lagrangian 

contains the term 

A^A(r) 
9 

which, when averaged over the cell, gets non vanishing contribution only from the diagonal terms in the double sum 
over the plane waves. 

The energy integration is performed by the residue theorem and the phase space is divided into different regions 
according to the pole positions, defined by 



, £ 2 + A |(v,e) . (40) 
Therefore we get 

PAln^ = f [ [ y A «( v ' e > -±[[ - kA (41) 

Ao J Pk 4vr v /e 2 + A^(v,e) tiJ J Pk ^ y/e + fc 2 A 2 

where the regions Pk are defined as follows 

P k = {(v,0\A E (v,e) = kA} (42) 

and we have made use of the equation 

(43) 

99 A v ' 

relating the BCS gap Ao to the four fermion coupling g and the density of states. The first term in the sum, 
corresponding to the region Pi, has P equal contributions with a dispersion rule equal to the Fulde and Fcrrel case. 
This can be interpreted as a contribution from P non interacting plane waves. In the other regions the different plane 
waves have an overlap. Since the definition of the regions Pk depends on the value of A, their determination is part 
of the problem of solving the gap equation. 

Let us comment on our result. We have shown that in our approximation the dispersion relation for the quasi- 
particles has several branches corresponding to the values fcA, k = 1, • • • , P. Therefore the following interpretation of 
the gap equation l|41l) can be given. Each term in the sum corresponds to one branch of the dispersion law, i.e. to 
the propagation of a gapped quasi-particle with gap fcA. The corresponding region is nothing but Pk- However, the 
regions Pk do not represent a partition of the phase space since it is possible to have at the same point quasi-particles 
with different gaps. 

As a final remark, it is possible to test our method in an exactly solvable case, i.e. the particular case of P identical 
plane waves. In this case the only non empty region in the gap equation is Pp therefore the gap equation reduces to 
the FF gap equation for a gap PA, i.e. the expected result. We assume this as a consistency test of our approximation. 

IV. NUMERICAL ANALYSIS 

The free energy Q is obtained by integrating in A the gap equation. At fixed Sfi, is a function of A and z q , 
therefore the energetically favored state satisfies the conditions 

and must be the absolute minimum. We now analyze three different structures by this criterion. We will assume 
[i = 400 MeV and S = /z/2. 
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A. Strip 

This is the case of two antipodal plane waves, i.e. we take P = 2 in Ea. (|35|l and 

n 1 = (0,0,+l), n 2 = (0,0,-1), 

in Eg. 1)36(1 . If e is the pole position in the £q complex plane, the gap equation is 



2A In ■ 



2d 



dv A B (v,e) 
-df 



r/v 



-de- 



A R (v,e) 



A o J J Pi ^ V^ 2 + A|(v,e) J Jp 2 4tt v /e 2 + A|(v,e) 
Here / represents a region where the two A e // appearing in A# have no overlap. In this region 
A E (v, e) = A + f^)^) - 20(2£)0(2£)0(i£)0(2!S)} A = A . 



(45) 



(46) 



In other words 



dv A_b(v, e) 



Pl 4vr ^/e 2 + A|(v,e) 



dv (g(^)g(^) + e{El)6{El) 2e{El)6(E\)6{El)6{El)} l 
— — dt 

4tt ^/e 2 + A 2 



(47) 



(48) 



On the other hand the integral over P 2 represents the region where the two A e f t appearing in A^ do overlap. In this 
region A|,(v, e) = 4A 2 and we have 



dv _ d £. A B (v,e) 



p 2 4^ vF+A|(^y 



dv 2{8(El)8(Ei)e(El)8(E2)} 2A 
4tt VC 2 + 4A 2 



(49) 



where the subscript 2A on the r.h.s. means that in the dispersion laws E^ 2 d one has to use Eqns. (J5J with A — ► 2A. 
As already observed, a test of our approximation is given by the fact that the gap equation obtained in this way 
reproduces in the case of qi = q 2 the FF result. 

As we observed above, we expect more reliable results for large A, i.e. far off the GL point. Nevertheless we can 
determine the phase transition point 6/i2 within this approximation; we find <5/i 2 ~ 0.83Ao, which compares fairly 
well with the exact result in the weak coupling limit, i.e. 0.75Ao; in our approximation we find a first order phase 
transition (see subsection IIV Cl for further comments). 

It may be useful to have a pictorial representation of the various regions. They are reported for the strip in a 
(z = cos#,£) plane in Fig. [2 at S/j, = 6/j,i, where our approximation is expected to work better. 




0.5 1 



FIG. 2: The regions Pi and P2 for the strip in the (z, £) plane with z = cos 8 — v ■ qi. The region Pi is the dark region, 
whereas the region P2 covers all the plane and is in light grey. Here A = 0.75Ao, z q = 1.0, <5/i = <fyn; £ in MeV. 

To get a better insight of the main features of our approximation, we can compare the approximate dispersion law 
we are using with the exact one. The exact dispersion law for the strip was derived in (see also [2{|; in both papers 
the energy is computed for each species from the corresponding Fermi surface): 
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El td = ±6nT<l-y + w\S\, El d = ±5 f ,±c l -v + w\^\ (50) 

where 

w- 1 = h (-^) ■ (51) 
\ q cos J 

Io is the Bessel function and 9 the angle between the z— axis and q = qi- 

This comparison for the strip is reported in Fig. [3](we use the parameters A = 0.75Ao, Sfi = Sfj,i and z q m 0.8, which 
is more appropriate for comparison with the exact result). On the right diagram the gray region depicts the area where 
coupling is more favorable, because here both quasiparticles with gap A and 2A contribute; to this area corresponds 
on the left diagram, in gray, the pairing region, where £* d > (fc = 1,2). The white areas depict respectively, on 
the right an area where coupling is less favorable, because only quasiparticles with gap 2A can contribute, and on the 
left the blocking region, where the energies are negative. We have limited the analysis to the region of smaller £ since 
this is the region more relevant for coupling. The two graphs do not coincide, which is due to the approximation, 
however they share some features, e.g. for both the central region (z ss 0) is disfavored; also the areas with small £ 
and large | cos(9| are less favorable in both graphs. This leads to the following interpretation of the different regions 
of the phase space. The region where only the quasiparticles of gap 2A are present corresponds grosso modo to the 
blocking region of the exact dispersion law, whereas the region where both species of quasiparticles (i.e. with gaps A 
and 2A) are present corresponds to the pairing region of the exact dispersion law. This interpretation is supported by 
the numerical results for the gap equation. In fact we find that the contribution of the former region (corresponding 
to blocking) to the gap equation is only 20% of the total. 




-1 -0.5 0.5 1 -1 -0.5 0.5 1 

z z 



FIG. 3: On the left: Pairing (in gray) and blocking (in white) regions as computed by the exact dispersion law for the 
strip. On the right: The gray region depicts the area where coupling is more favourable (here both quasiparticles with gap A 
and 2A contribute); in the white area only quasiparticles with gap 2A can be present. On the horizontal axis z = cos#; on the 
vertical axis the longitudinal momentum measured from the Fermi surface, £ (in MeV). Here A = 0.75Aq z q — 0.8, S/j, = Sjjbi. 



B. Octahedron and Face Centered Cube 

For the octahedron (i.e. body-centered-cube, bcc) we have P = 6 in (|35|l and we take 

ni = (+1,0,0), n 2 = (0,+l,0), n 3 = (0, 0, +1), 

n 4 = (-1,0,0), n 5 = (0,-1,0), n 6 = (0,0,-1) (52) 
in 136|) . For the face centered cube (fee) we have P = 8 and 

ni = -^(+l,+l,+l), na = i(+l,-l,+l), 
iib = 4t(-1,-1,+1), n 4 = 4e(-1, +!,+!), 
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n 5 = -J=(+l+,l 1 -l) 1 n 6 = -L(+l,-i,-i), 
n 7 = -^=(-1,-1,-1), n 8 = -^(-l,+l,-l) . 



(53) 



The topological structure of the different regions Pk is more involved for the cubes than for the strip. We simply offer 




FIG. 4: The pairing regions Pi, . . . , Pg (from left to right and from top to bottom) for the face centered cube in the (z, £) plane 
for a fixed value of the polar angle (p = .5 rad. In each figure the grey area corresponds to the pairing region. Here Sfi = S^ii, 
A = 0.21Ao, z q = 0.9. Note that for this particular value of the parameters the Pg region covers all the Fermi surface. 

in fig. 0]a picture of the different pairing regions Pi, . . . , Pg for the face centered cube in the (z,£) plane at fixed polar 
angle (p = 0.5 rad. In each figure the grey areas correspond to the pairing region. 



C. Numerical results 



We now present numerical results for four crystalline structures (FF, strip, bcc, fee) that are obtained solving 
Ea. (|38fl with the appropriate set of unit vectors n m . In Table [I] we report the results obtained at the Clogston point 
(Sfj, = S/ii). Here P is the number of plane waves of each structure; the table shows that, among the four considered 
cases, the favored structure is the octahedron [P = 6). 



p 


z q 


A 
A 


2sn 


1 


0.78 


0.24 


-1.8 x 10~ 3 


2 


1.0 


0.75 


-0.08 


6 


0.9 


0.28 


-0.11 


8 


0.9 


0.21 


-0.09 



TABLE I: The gap, z q = 8fi/q and the free energy at 8fi — <5/ii = Ao/v2 for different crystalline structures. 

We have analyzed also the case (5/i 7^ 6/ii. In fig. [5] we report the plot of the free energies of the octahedron (dashed 
line) and of the fee (full line) phases as a function of 5/i/Ap. 

We find that the octahedron is the favored structure up to 8^, « .95Ao. For larger values of Sfj, < 1.32Ao the favored 
structure is the fee. In Table ITT1 we report numerical results for 5/J.2 for each crystalline structure and the computed 
order of the phase transition between the LOFF and the normal phases. We have reported also the values of z q and 
of the discontinuity in A/Ao at = 6^2 — + . The values of z q and A are determined minimizing the free energy. 

As remarked above, for the strip, in our approximation neither the order of the transition nor the point where the 
transition occurs coincide with those obtained within the Ginzburg-Landau approximation [lj. The difference in 8^,2 
is ~ 10% and in z q is ~ 17% (the value for z q is 0.83). 
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FIG. 5: The values of the free energies of the octahedron (dashed line) and of the fee (full line) crystalline LOFF structures as 
a function of 5/i/Ao. The octahedron is the favored structure up to Sfj, « .95Ao; for .95Ao < 5/j, < 1.32Ao the fee is favored. 
Here, for each value of Sfi, the values of z q and A are those that minimize the free energy. 



p 


<5/te/Ao 


Order 


z q 


A/Ao 


1 


0.754 


II 


0.83 





2 


0.83 


I 


1.0 


0.81 


6 


1.22 


I 


0.95 


0.43 


8 


1.32 


I 


0.9 


0.35 



TABLE II: The values of <5/i2, z q — S/j,/q, the discontinuity of A/Ao and the order of the phase transition between the LOFF 
and the normal phases for different crystalline structures. 

V. DISCUSSION AND CONCLUSIONS 

In this paper we have discussed the gap equation for different Fulde-Ferrel-Larkin-Ovchinnikov phases at T = 0. 
Generally speaking this is a quite complicated problem and the only soluble case corresponds to the original Fuldc- 
Ferrel phase with a single plane wave. Typically the discussion of the gap equation is made via the Ginzburg-Landau 
app roximation which, however, holds only in proximity of a second order transition. However, the analysis of Ref. 
|14j shows that in many interesting cases there is no second order transition to the normal phase. Therefore, in this 
paper, we have developed a new kind of approximation to overcome the problems of the Ginzburg-Landau treatment. 
In fact, our approximation holds in the opposite case, that is when the gap is not too small. The main feature of our 
scheme is a convenient average over the sites defined by the crystalline structure of the condensate. The description 
one obtains amounts to the definition of a multi-valued gap function having P branches. Each of the branches has a 
gap given by fcA, k = 1, • • • , P, with A the solution of the gap equation. Furthermore each non vanishing value of the 
gap defines a particular region in the phase space, the pairing region Pk- It is interesting to notice that each region Pk 
depends on the actual value of A and therefore its very definition depends on the solution of the gap equation itself. 

We have noticed that the regions Pk do not represent a partition of the phase space since it is possible to have at 
the same point quasi-particles with different gaps. This observation would allow to formulate the problem in terms 
of an effective Lagrangian containing P velocity-dependent fields ip( k > , each of them defined in the region Pk, and 
therefore with gap fcA. We will come back to this effective Lagrangian in future work. 

The gap equation has been solved numerically. We have first analyzed the case P = 1 where our approximation is 
exact, and then the cases P = 2 (two antipodal plane waves, or strip), P — 6 (the octahedron or body-centered-cube) 
and P = 8 (the face-centered-cube) . Our main result is that, at T = 0, and up to <5/i « 0.95Ao the energetically 
favored crystalline structure is the octahedron (P = 6). For larger values of S/j,/ Aq the fee becomes the most favored 
structure up to the transition to the normal state occurring at 6[i2 = 1.32Ao- 

Also, in all the cases P = 2, 6, 8, we have found that the transition to the normal state is a first order one. 

It is difficult to compare our results to previous work since we have studied a region of the phase diagram not 
accessible to the Ginzburg-Landau approximation. In any event let us summmarize some of the previous results. 
Within the GL approximation, the authors of Ref. Q] found that the strip was favored in comparison to the one- 
plane-wave. Also in Ref. 0] a GL analysis is performed. These authors study a much larger number of crystalline 
structures and conjecture that the most stable case is the face-centered cube; however no second order transition is 
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observed for this structure and this makes questionable the Ginzburg-Landau approximation. As for the octahedron, 
the phase transition between the LOFF and the normal phase is also first-order. These results on the order of the 
transition for P — 6, 8 agree with our approximation. Finally, in a recent paper [22, the analogous problem has been 
studied in the context of the quasiclassical equations for superconductors by using a Fourier expansion. These authors 
get that the favored state at zero temperature is the octahedron, P = 6, and also find that the transition to the normal 
state is first order. Therefore the results of their analysis agree with our approximation, at least qualitatively. 
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